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1. Introduction 

Many eigenvalue correlations for random matrix ensembles as matrix Green functions, 
the /c-point correlation functions PQ [2] as well as the free energy [3] can be derived 
by generating functions. These functions are averages over ratios of characteristic 
polynomials for random matrices. 

A common approach to calculate generating functions is the supersymmetry method 
[H El HI E]- Another approach is the orthogonal polynomial method [7j. In the 
supersymmetry method one maps integrals over ordinary matrices to integrals over 
supermatrices. Its advantage is the drastic reduction of the number of integration 
variables. Nevertheless there is a disadvantage. Up to now it is not completely clear how 
to get the full structures found with the orthogonal polynomial method for factorizing 
probability densities. For such probability densities the /c-point correlation functions can 
be written as determinants and Pfaffians of certain kernels. This property was extended 
to the generating functions [H [9j [TUl E] • Unfortunately, the determinants and Pfaffians 
were not found for the full generating function after mapping the integrals over ordinary 
matrices to integrals over supermatrices. Only the determinantal expression of the k- 
point correlation function for rotation invariant Hermitian random matrix ensembles 
could be regained with help of the supersymmetry method [H [T2]. Gronqvist, Guhr 
and Kohler studied this problem from another point of view in Ref. [5]. They started 
from the determinantal and Pfaffian expressions of the /c-point correlation functions for 
Gaussian orthogonal, unitary and symplectic random matrix ensembles and showed that 
the kernels of the generating functions with two characteristic polynomials as integrals 
in superspace yield the known result. 

Changing coordinates in superspaces causes serious problems since the Berezinian 
rrS] playing the role of the Jacobian in superspace incorporates differential operators. 
These differential operators have no analog in ordinary space and are known in the 
mathematical literature as Rothstein's vector fields [14] . For supermanifolds with 
boundaries such differential operators yield boundary terms. In physics these boundary 
terms are called Efetov-Wegner terms [T5], [161 El HH]- They can be understood as 
corrections to the Berezinian without Rothstein's vector fields. There were several 
attempts to explicitly calculate these vector fields or the corresponding Efetov-Wegner 
terms for diagonalizations of Hermitian supermatrices but this was only successful 
for low dimensional supermatrices, e.g. (1 + 1) x (1 + 1) Hermitian supermatrices 
[19j [201 EI] • For higher dimensions the calculation of all Efetov-Wegner terms becomes 
cumbersome. Also for other sets of supermatrices such changes in the coordinates were 
studied [221 E31 EI] . There are quite general formulas [HI E21 [25] for calculating the 
Efetov-Wegner terms but for concrete examples these formulas become cumbersome and 
one never came beyond the low dimensional cases. As an example for such an approach 
we tried to derive all Efetov-Wegner terms by considering these terms as boundary 
terms resulting from partial integrations of differential operators which are equivalent 
to the integration over the Grassmann variables [18]. Though a quite compact form 
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for a differential operator for the Hermitian supermatrices was found we were unable 
to calculate all Efetov-Wegner terms. The order of such a differential operator is the 
number of pairs of Grassmann variables. 

We derive the supermatrix Bessel function [25] for Hermitian supermatrices with 
all Efetov-Wegner terms. Thereby we apply a method called "supersymmetry without 
supersymmetry" [10] on Hermitian matrix ensembles whose characteristic function 
factorizes in the eigen-representation of the random matrix. This approach uses 
determinantal structures of Berezinians without mapping into superspace. We combine 
this method with the supersymmetry method described in Refs. [H [T2] . The supermatrix 
Bessel function is then obtained by simply identifying the left hand side with the right 
hand side of the resulting matrix integrals. As a simple application we calculate the 
generating function of arbitrary, finite dimensional, rotation invariant Hermitian matrix 
ensembles with and without an external field. Hence we generalize known results 
[H [271 E21 ED] because they now contain all Efetov-Wegner terms which guarantee the 
correct normalizations. 

We organize the article as follows. In Sec. [2J we will give an outline of our 
approach and introduce some basic quantities. Using the method "supersymmetry 
without supersymmetry" we derive a determinantal structure for generating functions 
of rotation invariant Hermitian random matrices without an external field and with a 
factorizing characteristic function, in Sec. [3j In Sec. HI we briefly present the results 
of the supersymmetry method. We also calculate the supermatrix representation of 
the generating function with one determinant in the numerator as in the denominator. 
For this generating function the corresponding supermatrix Bessel function with all 
Efetov-Wegner terms is known [28] since it only depends on (1 + 1) x (1 + 1) Hermitian 
supermatrices. In Sec. [5j we plug the result of Sec. 0] into the result of Sec. [3] and 
obtain the supermatrix Bessel function with all Efetov-Wegner terms for arbitrary 
dimensional Hermitian supermatrices. Moreover we discuss this result with respect to 
double Fourier transformations and, thus, Dirac-distributions in superspace. In Sec. El 
we, first, apply our result on arbitrary, rotation invariant Hermitian random matrices 
without an external field and, then, in the presence of an external field. Details of the 
calculations are given in the appendices. 

2. Outline 

We want to generalize the well known result for Gaussian superfunctions [281 EE] 

f ( 1 G , 2 , G+ A ,r I , («1 ~ «g) f exp {-[sj + s|/2 + ISiKl + S 2 K 2 ) M 

J exp (--Stra + zStr™j d[a] oc 1 - J — d[ 

( 

The matrix a is a Hermitian (1 + 1) x (1 + 1) supermatrix with the bosonic eigenvalue 
si and the fermionic eigenvalue is 2 and k = diag (/Cj, K2) is a diagonal (1 + 1) x (1 + 1) 
supermatrix. The first term on the right hand side of Eq. ([T]) is the Efetov-Wegner term 
and the second one is the supergroup integral appearing after diagonalizing a = UsW 
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with U G U(l/1). Indeed, this result can be easily generalized to higher dimensional 
supermatrices with Gaussian weights. However, considering Gaussian weights is not 
sufficient to generalize this result to arbitrary weights. Thus, we pursue another 
approach. 

The main idea of our approach is the comparison of results for generating functions, 



p(N) {H 



k 2 

n det(H 



K j2 l N ) 



-d[H] , 



(2) 



Herm (N) 



II det(H 
i=i 



jl^N) 



obtained with and without supersymmetry. The matrix Ujv is the N dimensional unit 
matrix. The integration domain in Eq. (j2J) is the set Herm (N) of N x N Hermitian 
matrices and the parameter k = diag («n, . . . , k^i, «i2, • • • , Kktf) — diag (ki, k 2 ) are 
chosen in such a way that the integral exists. One common choice of them is as real 



energies Xj with a small imaginary part is and source variables J,-, i.e. k 



[TJ [12]. The differentiation with respect to Jj of zj^l^ generates the matrix Green 
function which are intimately related to the fc-point correlation functions [28]. The 
measure d[H] is defined as 

N N 



J 1/2 



x .i "F Jj 



-IE 



d[H] = J] dH nn H dRe H mn dlm H n 



(3) 

n=l l<m<n<N 

The matrix Ijy is the N x N unit matrix and P^' is a probability density over N x N 
Hermitian matrices. 

In Sec. El we show that Z^V k2 is a determinant, 



4f/L(«)~det 



if the Fourier-transform 
FP {N \H) 



Z\J\ (l^alj K b2) 



h 



al 



K b2 



P( N \H) exp 



itrHH 



d[H] 



(4) 



(5) 



Herm (TV) 



factorizes in the eigen-representation of the Hermitian matrix H. These determinantal 
structures are similar to those found for generating functions with factorizing p( N > in 
the eigen-representation of the Hermitian matrix H [91 [10] . 

Moreover, the integral of Z^l can be easily mapped into superspace if the 
characteristic function JFp( N > is rotation invariant, see Refs. [U [12]. In this 
representation one does not integrate over ordinary matrices but over Wick-rotated 
Hermitian supermatrices [TS1[T2]. The integrand is almost rotation invariant apart from 
an exponential term which reflects a Fourier-transformation in superspace. We aim at 
the full integrand with all Efetov-Wegner terms appearing by a diagonalization of a 
supermatrix. In particular we want to find the distribution (pk 1 /k 2 which satisfies 



exp [-zStr Kp) F(p)d[p] = I Bexf\ k (r)cp kl / k2 (-ir,K)F(r)d[r 



k 1 /k 2 



(6) 
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for an arbitrary rotation invariant superfunction F. In Eq. ([6]) we diagonalize the 
(ki + k 2 ) x {ki + /c 2 ) Hermitian supermatrix p to its eigenvalues r. This diagonalization 
does not only yield the Berezinian Ber^^ but also the distribution tpki/ks- The 
distribution (fk 1 /k 2 is the integral over the supergroup V (ki/k^) and, additionally all 
Efetov-Wegner terms. It is called the supermatrix Bessel function with Efetov-Wegner 
terms [261 HE]. 

We derive <fki/k 2 i n two steps. In the first step we combine the mapping into 
superspace with the determinant (jlj) for factorizing FP( N \ This is a particular case of 
the identity Eq. but it is sufficient to generalize its result to an arbitrary rotation 
invariant superfunction F in the second step. 

The procedure described above incorporates determinants derived in Ref . [221 OH] ■ 
Without loss of generality let p > q. Then, these determinants are 

1] («al - «6l) II ( K a2 ~ K b2 ) 
l<a<b<p l<a<b<q 



Ber«(«) 



p i 

n n («ai 

a=l 6=1 



(7) 



«62 



,_ 1)P (p-l)/2 ^ 



^p-q^q-p 

Hi K a2 
K bl — K a2 



l<a<q 
l<b<p 



<P-Q 
l<b<p 



(8) 



-l) p(p - 1)/2 det 



1 



«61 — K a2 J i^ ^? 

l<6<p 

Kr 1 } 1<< 

l<a<p— q 
l<b<p 



All three expressions find their applications in our discussion. For p 
Cauchy-determinant [28] 

1 



(9) 



q we obtain the 



Ber; 2 / ;(«) = (-l)^- 1 )/ 2 det 



l<a,b<p 



whereas for q = we have the Vandermonde-determinant 

' (-l) p(p - 1)/2 det[^} 1<a 



(2) / 
p/0\ 



Ber , n (K) 



Ap(«i) 



(10) 



111 



These determinants appear as square roots of Berezinians by diagonalizing Hermitian 
supermatrices. This also explains our notation for them. The upper index 2 results from 
the Dyson index /3 which is two throughout this work. Thus we are consistent with our 
notation used in other articles [121 [30], [TOj [11] . 



3. The determinantal structure and the characteristic function 



We consider the generating function 02]). Let the probability density be rotation 
invariant. Then, the characteristic function ([5|) inherits this symmetry. We consider 
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such characteristic functions which factorize in the eigenvalues E = diag (Ei, . . . , En) 
of the matrix H, i.e. 

JV 

FP {N) (E) = Hf(E J ). (12) 
i=i 

Due to the normalization of P^ N \ the function / : R — > C is unity at zero. We restrict 
us to the case 

k 2 <k x <N. (13) 
This case is sufficient for our purpose. 



Let / be conveniently integrable and d — JV + k 2 — k\. We derive in Appendix A 
that the generating function in Eq. has for factorizing characteristic function (1121) 
under the condition ([TBI the determinantal expression 



7 (AT) , N (_l)fc2(fe 2 + l)/2+fc 2 fci 



Z^M k = v ^det 

.(2) 



'Ber2> M 



K al — | l<a<fci L ' J l<a<fci 

l<b<k 2 d+l<b<N 



(14) 

for all G {d, d + 1,...,A/"}. Notice that Eq. (1141 is indeed true for all values 
A^ G {d,d+ 1,...,N} because the determinant is skew symmetric. Hence, we may 
add any linear combination of the last k% — k 2 columns to the first k 2 columns. 

With this theorem we made the first step to obtain all Efetov-Wegner terms for 
the diagonalization of a supermatrix. For non-normalized probability densities, we find 
Eq. flH| multiplied by / A (0) where A = d(d + l)/2 - (JV - l)JV/2 - Nk 2 . We need the 
non-normalized version to analyze the integrals in Sec. El 

Equation (fl4"|) is also a very handy intermediate result. It yields a result of a simpler 
structure than the one for factorizing probability densities [9j [10]. Thus, all eigenvalue 
correlations for probability densities stemming from a characteristic function with the 
property f jT2|) are determined by one and two point averages. 



4. Integral representation in the superspace of the generating function and 

(N 
1/1 

In Refs. [U El [121 EQ] , it was shown that the generating function ([2]) can be mapped to 

)) 

i/k 2 

,t 



the Efetov-Wegner term of Z^V 



an integral over supermatrices. Let be the set of supermatrices with the form 



P 



(15) 



Pi V ] 

V P2 

The Boson-Boson block p\ is an ordinary k\ x k\ positive definite Hermitian matrix and 
the Fermion-Fermion block p 2 is an ordinary k 2 x k 2 Hermitian matrix. The off-diagonal 
block rj comprises k\ x k 2 independent Grassmann variables. We recall that (f]'y = —r\ 
and the integration over one Grassmann variable is defined by 

J lfnm d Vnm = J VnL dr ln m = fyl > J ^ {0, 1} . (16) 
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We need the Wick-rotated set Sj^L = n^E^y^n^ to regularize the integrals below. 
The matrix IT^ = (l kl , e^ /2 t k2 ) with i[) e]0, ?r[ is the generalized Wick-rotation [ISJ[T2]. 
We assume that a Wick-rotation exists such that the characteristic function is a 
Schwartz function on the Wick-rotated real axis. Examples of such functions are given 
in Ref. [12]. 

We define the supersymmetric extension $ of the characteristic function TP with 
help of a representation 

TP l (tTH m ,m eN) = TP {N \H) (17) 

function in matrix invariants, 

$(fri/feO( p ) =J^P 1 (Strp m ,mGN). (18) 

For a factorizing characteristic function (TT2]) we want to consider only such extensions 
<3>( fc i/ fc 2) which factorize, too. In particular, we only use the extension 

a=l 6=1 ^ v bl > 

Indeed the extension of TP^ to $( fc i/ fc 2) is not unique and, hence, the factorization 
property (TT9T) does not hold for each extension, see Ref. [12J. However, for our purpose 
the choice (119]) is sufficient at the moment. Later on we extend our result to arbitrary 
rotation invariant superfunctions such that it also applies to other rotation invariant 
extensions. 

Let the signs of the imaginary parts for all Kj\ be negative. Assuming that $ 
is analytic in the Fermion-Fermion block p%, the generalized Hubbard- Stratonovich 
transformation [TJ [12j [30] tells us that the integral Eq. (J2J) with the condition ([13]) is 

k2 / o \ d-1 



= C[ N X j $^(p) expHStr K p]det d Pl f[ (e"*^) (r, 2 ) d[p] , 

(20) 



k l/ k 2 



where e l ^Vj2 are the eigenvalues of P2 and the matrix p is given by 



P 



The constant is 



Pi 






e^(p2 + ??prV) 



(21) 



(_ 1 )fc 2 (fc 2 + 2 jv-i)/ 22 jv(fc 2 -fc 1 ) 7r jV(fc 2 -fc 1 ) + fc 22 fc 2fcl+ fc 2 -fc 1 vd ( U (jV)) 
fcl/fc2 ~ [(d - l)!] fe2 Vol(U(d)) / (0) ' (22) 

Here, the volume of the unitary group U (N) is 

2tt j ' 



Vol(U(JV)) = n 7 r- rf? . (23) 
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The term /(O) becomes unity if we consider normalized probability densities. The 
definition of the measure d[p] = d[pi]d[p2]d[r]} is equal to the one in Ref. [30J, 

fci 

d[pi] = 17 dp nn i ] J dRe p mnl dlm p mnl , (24) 

n=l l<m<n<ki 

d[p2\ = e tk ^ Yl dp nn2 Yl dRe Pmn2dlm p mn2 , (25) 

n=l l<m<n<k2 
fel k 2 

d[rj] =e-^Y[Y[d Vmn d V * mn . (26) 

n=l m=l 

We use the conventional notation for the supertrace "Str" and superdeterminant 
"Sdet ". 

Let p G Yy^\ k ^ with the form ( Fl5|) . As long as the eigenvalues of the Boson-Boson 
block pi are pairwise different with those of the Fermion-Fermion block p 2 , we may 
diagonalize the whole supermatrix p by an element U G XJ(ki/k 2 ). The corresponding 
diagonal eigenvalue matrix is r = diag (r 11; . . . , r fcll , e"V 12 , . . . , e"V 22 ) = diag (r 1: e^r 2 ), 
i.e. p = UrW. Due to Rothstein's [14] vector field resulting from such a change of 
coordinates in the Berezin measure, we have 

d[p] ^Ber ^ /k2 (r)d[r}dp(U), (27) 

where d[r] is the product of all eigenvalue differentials and dp{U) is the supersymmetric 
Haar-measure of the unitary supergroup U (ki/k 2 ). We have to consider some boundary 
terms since the Berezin integral is fundamentally connected with differential operators 

P31SDHBIE21. 

We consider the generating function zffl. Equation (TSUI) is an integral over Dirac 
distributions. Hence, we cannot simply apply a Cauchy-like theorem [161 HHJ E31 EH El 
[18] but also for this integration domain we obtain an Efetov-Wegner term for Z^V. Let 
the function 1/f be analytic at zero. Then, we have 



f N (0) 



exp[— tStr Kr] 



z(-l) 



N 



+ (N-1)\ 



r=0 

JV-1 



J(e^r 2 ) 

[ Kl - K2 J} Tl) r exp[-zStr *r]rf (e~^^-] 5 (r 2 ) dr 2 d n (2* 



-1) 



N 



(N-l)\ 



1 ( 9 +e ^d\ +i _K 1 -K 2 



r\ — e % ^r 2 \dr 1 dr 2 J T\ — e l ^r 2 



exp[-*Str «r]rf e"^— 5 (r 2 ) dr a dn . (29) 
In Eq. f[2"9"j) . we, first, integrate over r 2 and, then, over r x . We derive this result in 



N-l 



Appendix B The first summand on the right hand side of equality (1281) is 1 which is 



the Efetov-Wegner term. The second equality (1291) is more convenient than equality 
for the discussions in the ensuing section. 
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5. Supermatrix Bessel function with all Efetov-Wegner terms 



Let supermatrices in Yr^\ k2 be similar to those in Tr^\ h ^ without the positive definiteness 
of the Boson-Boson block. We want to find the distribution (f>k 1 /k 2 which satisfies 



exp [— zStr up] F(p)d[p] 



Ber 



kx/k 2 ( r )Vk 1 /k 2 (-w, K)F(r)d[r] 



(30) 



g(V0 M. k i+ k 2 

k 1 /k 2 

for an arbitrary sufficiently integrable, rotation invariant superfunction F analytic 
at zero. Recognizing that the integral expression f[2TJj) includes the supersymmetric 
Ingham-Siegel integral [TJ [T2J, [30] , the generating function is apart from a shift in the 
Fermion-Fermion block and analyticity a particular example of this type of integral. 

First we want to derive <fiki/k 2 in Eq. (1301) for factorizing superfunctions (JT9"j) and, 
then, extend to arbitrary F. We show in Appendix C that the distribution (fki/k2 defined 
by 



: «0 



F{p) exp[— «Str Kp]det d pi Y\ ( e 



■iil> 







dr 
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e-^5(r j2 )d\p) 



Be 4?/fc 2 ( r )WteH' 



k-2 



R^ 1 xR fc 2 

is given by 

(p kl/k2 (-ir, k) 



, K)F{r)dei d r x JJ ( e 

3=1 







dr 



J 2 



d-1 



(r j2 ) d[r] (31) 



(_l)(fcl+fc2)(fcl+fc 2 -l)/2^ 7r -)(/e2-fci) 2 /2-(fei+fc2)/2 



2^^/, 1 !A; 2 ! A /Ber( 2 i ) A2 ( K )Ber fci/fc2 



,(2) 



(r) 



x 5^ exp 

^2e6(fc 2 ) 



fci fc 2 



5=1 



(32) 



x det 
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.a-1 



+ e 



5 



<9r 



a; 2 (a)2 



l<a<fc 2 
K6<fci 



wi(6)l 



l<a<fci— fc2 
K6<fa 



for an arbitrary sufficiently integrable, rotation invariant superfunction F analytic at 
zero and factorizing in the eigen representation of p. The set of permutations over k 
elements is &(k). The derivation of Eq. (132]) incorporates Eqs. ()14p and ( |29j) . 

to arbitrary rotation 



In Appendix D we make the next step and generalize Eq. 
invariant superfunctions without the additional Dirac-distributions in the integral, cf. 
Eq. (|20|) . The distribution (fk 1 /k 2 defined by Eq. ( 130]) is indeed the one defined in 
Eq. ([32]) for an arbitrary sufficiently integrable, rotation invariant superfunction F(p) 
analytic at zero. For such superfunctions and for generalized Wick-rotations ip g]0, 7r[ 
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^* 1 /fca( -zr > K ) 



.]_)(fci+A 2 )(fci+fc2-l)/2j 



Ak 2 -k 1 ) 2 /2-{k 1 +k 2 )/2 



2 fclfc ^i!^!A/Berg ) /fc2 ( K )Ber^ /fc2( r) 



(2) 
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Wl£S(fci) 
l^2€6(fc 2 ) 



x det 



-2vre # 5(r Wl(fe)1 )5(r W2(a)2 ) exp (-z/Cbi^^i + ^2^(0)2) , 

H — XV^&l ~~ K a2j 



K&l - «a2 ^(6)1 - 6^^(0)2 



<a<fci— fc 2 
l<b<fei 



with the distribution 



, x = 0, 

1 , else. 



(34) 



Equation fl33|) is true because of the Cauchy-like theorem for (1 + 1) x (1 + 1) Hermitian 
supermatrices, see Refs. [TBI [3U [T3]. One has to careful on which half of the complex 
plane the general Wick-rotation is lying. If ip e]7r, 2ir[ then the minus sign changes to a 
plus infront of the Dirac-distributions. 

We emphasize that the distribution (pki/k 2 ( r i K ) is n °t symmetric in exchanging its 
arguments r and k. Apart from the characteristic function x such a symmetry exists 
for the supermatrix Bessel functions j2BHsS3HB] which is (pk 1 /k 2 ( r , K ) without the Dirac- 
distributions, i.e. 



-zr, k 



(_l)(fci+fc2)(fci+fe2-l)/2^ 7r ^(fc 2 -fci) 2 /2-(fci+A:2)/2 



X 



E det 



wies(fci) 
w 2 e6(fc 2 ) 



exp (-««6i»" Wl (6)i + «e # K a2 r W2(a) 2) 



2^^fc 1 !fc 2 !JBer^ /fc2 («)Ber fc7/fc2 (r 



(35) 



X(«61 - ^ a2 ) 



r Ul (b)i - e^r W2(a)2 

C(l)i eX P (-^1^(6)!)}^^ 

(_l)fc2(fc2-l)/2+fcifc 2 ^ 7r ^(fe 2 -fci) 2 /2-(fei+fc 2 )/2 



l<a<k 2 
l<6<fci 



-fca 



K6<fa 



n 



K al — Kb2j 



det [expHKairti)]^^*, det [exp(ze^« a2 r 62 )] x < ai6 < fe2 



l<a<fci 
l<b<k 2 



Ber Sfe( r ) 



The asymmetry, ^/^(r, k) 7^ (pk 1 /k 2 ( K , r ), is mainly due to the diagonalization of p 
to r whereas the supermatrix k is already diagonal. The characteristic function x m 
Eqs. f )33|) and f l35]) is crucial because of the commutator 

d . _„,, d 



dr 



+ e 



-iip 



wi(6)l 



exp (-iKwr^fm + «e^K a ar Wa ( a ) 2 ) 



*(«61 - «o2) ex P (-*«ML»«i(&)l + ie* K aa r Wa ( o ) 2 ) x(«u. - K a2 ) . (36) 
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Indeed the set which is cutted out by x is a set of measure zero and does not play any role 
when one integrates tp kl /k 2 {— ir i K ) or ¥k 1 /k 2 {~ ir i K ) over K with a function continuous 
around this set. However it becomes important for the integral 

Wfei/ka r)(p kl/ka (-ir, £)Ber [ 2) i/k2 (r)d[r\ (37) 

l<6<fc 2 V ' V / 

where s = diag (s u , . . . , s kll , e~^s m . . . , e~^s fa2 ) and k = diag (/t u , . . . , K kll , 
e~ l ^Ki2, • • • , e~^Kk 2 2) with s a b, K a 6 e R. This result is the correct one for the supermatrix 
Bessel function. The difference to other results ESI EES] is the distribution \ 
which guarantees that the Dirac-distribution (157)) in the eigenvalues s of a Hermitian 
supermatrix cr vanishes if a bosonic eigenvalue of k equals to a fermionic one. On the first 
sight this would also happen if we omit x due to the prefactor J Ber f^ k2 (k) in Eq. fl3Tl) 
which vanishes on this set. Yet after applying the distribution (157)) on functions over 

(2) 

the set of diagonalized Hermitian supermatrices comprising the Berezinian Ber k \ 2 (s) 
in the measure this term cancels out, i.e. 

F(k)= / ~ ^-Ber I > /k2 (s)d[s] 

^F(k) J] - e-^Ksa) ■ (38) 

l<a<fci 
l<b<fc 2 

The reasoning why Eq. (157)) has to vanish on this set becomes clear when we interpret 
Eq. (|57|) as an integral over the supergroup U(fci//c 2 ), i>e. 

^i/fc 2 (^> r )^fci/fc2(-^^) Ber i 2 1 ) /fc 2 ( r ) rf H ~ / S(UsW -K)dfx(U) . (39) 



fcl/fe2 ' 

K fe i+ fc 2 U(fei/& 2 ) 

The measure dfi(U) is the Haar measure on U(/ci//c 2 ) and the Dirac-distribution is 
defined by two Fourier transformations 

5{UsU ] - k) ~ J exp [iStr p(tf sC/ f - k)] . (40) 

The Haar measure d// of the supergroup U(fci/fc 2 ) cannot be normalized as it can be 
done for the ordinary unitary groups since the volume of U (/ci//c 2 ) is zero for k\k 2 ^ 0. 
This is also the reason why Eq. ( 157)) has to vanish if one bosonic eigenvalue of k equals 
to a fermionic one. Then the integral ( 139)) is rotation invariant under the subgroup 
U(l/1) which has zero volume, too. This cannot be achieved without the distribution 
X as it was done in the common literature [25J |55J Q2]. Interestingly the replacement of 
Vk 1 /k 2 by y2fci/fc 2 m Eq. ( 157)) yields the full Dirac-distribution 

£fci/fc 2 {is, r)(p kl/k2 (-ir, £)Ber {r)d[r} (41) 
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_ Trfe-fci) 2 det [5(K al - s bl )}^ a ^ ki det [e^8(n a2 - Sft 2 )] 1 < a|ft < fc2 



which is similar to the one in Eq. fl37|) . Equation ( 14 ip is derived in Appendix E 

We want to finish this section by a remark about the relation of the result for the 
supermatrix Bessel function, see Eq. ( |33|) . and the differential operator derived by the 
author in an earlier work [T5]. This differential operator is defined by 

D^F(r) = [ F(p)d[rj] (42) 



for an arbitrary sufficiently integrable superfunction F on the {k\ + k 2 ) x {k\ + k 2 ) 
supermatrices invariant under U {ki/k 2 )- It has the form 

D {klk2) = 7 (43) 

(fc 1 fc 2 )!(47r)^^A fcl (r 1 )A fcl (e^r 2 ) 1 ' 



n=0 \ J x / l<a<fci / 

l<6<fc 2 



2 



where we define 



5 2 fa 9 2 



str 9r2 - E 5r 2 e " 2 ^ E 9r 2 • ( 44 ) 

0=1 al 6=1 62 

Due to Eq. (142 p the differential operator Dr 1 ^ is equivalent to the integration over all 



Grassmann variables of the supermatrix p. 

The comparison of Eqs. ( H2]) and (j4"3~P with Eqs. (T30T) and (T33|) for an arbitrary 
sufficiently integrable, rotation invariant superfunction F and arbitrary diagonal 
supermatrix k yields 



(45) 



n(fafa) = 1 * TT f JL , P -t^ JM , /r pt - ( 2 ) fr^l 

r (27r)fefeA fex (ri)A fex ( e ^r 2 ) 1 Al l^ai dr b J V fo/feA^- 

X /C^ 

l<b<fc 2 

Thus, we have found a quite compact form for D^ lk2 ^ which is easier to deal with as 
the one in Eq. 



6. Some applications for Hermitian matrix ensembles 

In random matrix theory generating functions as 

Herm(AT) J_1 

are paramount important since they model Hermitian random matrices in external 
potentials [361 EI] or intermediate random matrix ensembles [371 EH EH1 HQl HD S2]. 
The matrix H is & N x N Hermitian matrix and can be arbitrarily chosen or can also 
be drawn from another random matrix ensemble. The external parameter a G K. is 
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the coupling constant between the two matrices H and Hq and yields the generating 
function (j2j) for a = 0, i.e. Z^ n \k,0) = Z^(k). The variables Kj% have to have an 
imaginary increment to guarantee the convergence of the integral, i.e. Kji = Xj% — Jj — ie. 

In subsection 16.11 we consider the generating function (|46| with A = 0. We will 
use the mapping of this integral to a representation in superspace which is shown in 
a previous work by the authors [10] and diagonalize the supermatrix. In a formalism 
similar to the case a = we will treat the more general case a ^ in subsection 16.21 

6.1. Hermitian matrix ensembles without an external field (a = 0) 

Omitting the index N in p( N ^ we consider a normalized rotation invariant probability 
density P with respect to N x N Hermitian random matrices. 

Following the derivations made in Refs. [HH2] we have to calculate the characteristic 
function J-P, see Eq. (jSJ). Assuming that this can be done we recall the rotation 
invariance of J-P. This allows us to choose a representation of TP as a function in a finite 
number of matrix invariants, i.e. J r P(H) = J r Pi(trH, . . . , tr H ). A straightforward 
supersymmetric extension $ of the characteristic function is the one in Eq. (fTT|) and its 
Fourier transform is 

J^$(a) = 2 2fc(fc - 1} J $(p)exp(-zStrp(7)d[p]. (47) 

The Fourier-transform is denoted by Q in Ref. [lj. The supermatrices p and a are 
Wick-rotated with the phases e 1 ^ and e -4 ^, respectively. 
The supersymmetric integral for zl N \K,aH ) is [TJ | 



Z ( k N) {K,aH )= J .F$(<7)Sdet ~V ® lj\r + al*+* ® #o -«)<*[*] ■ (48) 



y(-V-) 
^fc/fc 



Setting the coupling constant a to zero we also find 

z[ N \K)=2 2k ^ j J ^(a)ex P [tStip(a-K)]I N (p)d[a]d[p], (49) 

v«0 y,(-i>) 
^k/k ^k/k 

where the supersymmetric Ingham-Siegel integral is [lj 
I ( k N \p) = J Sdet~ N (a + iet k+k )exp(-iStr pa + eStrp)d[a} 



?(-</>) 
^k/k 



k 

^_ 1 ^k(k+l)/2 2 -k(k-l) JJ 
3=1 



N-1 



(50) 



The distribution B is the Heavyside distribution and Tji and e^r^ are the bosonic and 



fermionic eigenvalues of p, respectively. In Appendix F we perform the integration (I50p 
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with help of the results in Sec. [5] and find 



2fc(fc-i) 



d[s) 



(51) 



x det 



e^5(s bl )8(s a2 ) ( K a2 \ N N (e ^s a2 - K a2 ) 1 x(« 6 i - Waa) 



+ 



«w - K a2 V k w y ' 2vn (e-^s w - n bl ) N+1 (s bl - e~^s a2 



l<a,6<fc 

The first term in the determinant is the Efetov-Wegner term whereas the second term 
can be understood as integrals over supergoups. 



6.2. Hermitian matrix ensembles in the presence of an external field (a ^ 0) 

For Hermitian matrix ensembles in an external field it is convenient to consider the 
integral representation 



Z { k N) (K,aH ) = 2 2k{k - 1) J $(p)exp[-zStr«;p] 



(52) 



^k/k 



J exp(— «Str pa + £:Str p)Sdet 1 (cr <8> 1 N + atk+k ® H + ze:IL A r (fc+fc))'^[c r ] 



f (-0) 
_ k/k 



dip] 



for the generating function, see Eqs. ( 14~6]) and (T48|) . In Appendix G we integrate this 
representation in two steps and get 

(_l)fc(fc-i)/2 



zl N) (K,aH ) 



A N (aE ) x /Ber$ k (K) 



(53) 



det 



R 2k 



{B l (r b i,r a2 , n bl , n a2 )} l<a b<k {B b2 (r a2 , K a2 )} 



l<a<k 
Kb<N 



b<N 



Kb<k 



<&(r)d[r] 



where 



Bi(r bll r a2 ,K bl ,K a2 ) 



e-^5(r bl )5(r a2 ) 



+ % 



K>b\ — K a2 

l[ N \r m e^r a2 ) 



2n(r bl - e^r a2 ) 



exp (-iK bX r bl + ie l ^K a2 r a2 ) - K a2 ) (54) 



B b2 (r a2 , n a2 ) = exp (ze^K a2 r a2 ) ( -ze~ # -^— J e" iV, 5(r a2 ) JJ - K a2 ) (55) 



B 3 (r bl ,K bl ,aE^) 



- i exp (-^K&ir 6 i 



(laE^r^" 



71= N 



nl 



Y[ x ( Kbl - e ^K j2 ) .(56) 
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The case a = can be easily deduced from the result fl53|) 




x exp (— tStr nr) $(r)<i[r] . 



ir } (r bl ,e^r a2 ) 
27r(r w - e^r a2 ) 



J l<a,6<fc 



(57) 



Again we are able to distinguish the Efetov-Wegner terms from those terms 
corresponding to supergroup integrals. In the determinant of Eq. ( |57|) and in the left 
upper block of Eq. (153|) . see Eq. (154"]) . the Dirac-distributions are the contributions from 
the Efetov-Wegner terms. When we expand the determinants in the Dirac-distributions 
we obtain the leading terms of Z^ N \ zj^_\, . . ., which are exactly those found in 
Refs. [HH2]. Thus, we have found an expression which can be understood as a generator 
for all generating functions Z^ N \ 

The external matrix H can also be drawn from another random matrix ensemble as 
it was done in Refs. [371 EHl USE EH El E2] • However, we do not perform the calculation, 
here, since it is straightforward to those in an application of Ref. [TO] . 

7. Remarks and conclusions 

We derived the supermatrix Bessel function with all Efetov-Wegner terms for Hermitian 
supermatrices of arbitrary dimensions. We arrived at an expression from which one can 
easily deduce what the Efetov-Wegner terms are and which terms result from supergroup 
integrals. With this result we showed that the completeness and orthogonality relation 
for the supermatrix Bessel function without Efetov-Wegner terms slightly differs from 
the formerly assumed one [38J. It has to be zero on a set of measure zero and, thus, 
does not matter for smooth integrands but it plays an important role if the integrand 
has singularities on this set. 

We applied the supermatrix Bessel function with Efetov-Wegner terms to arbitrary, 
rotation invariant Hermitian random matrix ensembles with and without an external 
field. The already known leading terms [33 ESI EHl ESI SOI IH1 HI [23 S21 H21 HD] were 
obtained plus all Efetov-Wegner terms. The correction terms were unknown before in 
this explicit form and yield new insights in the supersymmetric representation of the 
generating functions. In particular the Efetov-Wegner terms become important for the 
matrix Green functions. 

We also found an integral identity for the generating functions whose integrand 
can be easily expanded in the Efetov-Wegner terms. In such an expansion one obtains 
correlation functions related to fc-point correlation functions which are of lower order 
than those corresponding to the originally considered generating function. Thus, it 
reflects the relation of Mehta's definition [7] for the fc-point correlation function and the 
one commonly used in the supersymmetry method [JJ which was explained in Ref. [43J . 

We expect that similar results may also be derived for other supermatrices, e.g. 
diagonalization of complex supermatrices [44] . Nevertheless we guess that the knowledge 
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about the supergroup integrals as well as about the ordinary group integrals is crucial. 
We could only obtain these compact results due to this knowledge. 
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Appendix A. Derivation of Eq. (114ft 

We plug the characteristic function fl5]) in Eq. (j2J) and diagonalize H. This yields 



k 2 



7 (N) , x _ Vn 
^fcl/faW ~~ 2 N 7T N ' 2 



exp 



-ttr HE 



H det(H - K j2 l 



N) N 



Herm (AT) R N 



where the constant is 



N 



Vi 



N 



N\ 11 (j - 1)! 



II m)A%(E)d[E]d[H) , 

fl det(H - Kjil N ) i= l 

(A.1) 



(A.2) 



The diagonalization of H yields the matrix Bessel function [15J H6] according to the 
unitary group U (JV), 



Vn(E,E) 



exp 



:tl EUElP 



d(j,(U) 



V(N) 
N 



det 



exp(-tE a E b ) 



l<a,b<N 



fi A N (E)A N (E) 
The measure d/i is the normalized Haar-measure. Thus, we find 



(A.3) 



Z (N) (k) 



t N(N-l)/2 
~N~ 



det 



exp(-iE a E b 



Ka,b<N 



3=0 

k 2 

N EI (Ea — K b2) 

x Hf(E a ) b -^ A N (E)A N (E)d[E]d[E\. 



(A.4) 



a=l 



U(E a -K bl ) 
6=1 



Please notice that we integrate first over the variables E and then over E. Here, one 
can easily check that the normalization is Z^ k ^(0) = f N (0). Since determinants are 
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skew-symmetric, we first expand the Vandermonde determinant An(E) and then the 
determinant of the exponential functions. We have 



N 



(2*0*110' -1)! 



N 

n 

a=l 



f(E a ) exp(-zE a E a )E a a - lb -^ 



n (#«-««) 

6=1 



x A N (E)d[E]d[E] . (A.5) 

Following the ideas in Ref. [TO], we extend the integrand by a square root Berezinian 
and find with help of Eq. ([9]) the determinant 

(_l^k 2 (k 2 ~l)/2+(k 2 +l)k H N{N-l)/2 i 



Z {N) (K 



jV 



(27r)*n(j'-l)! 



x det 



l<6<fc 2 



l<a<d 
l<b<k 2 



2 

f{E 2 )E b 2 ' 1 E a 1 - 1 exp[-iE 2 E 1 \d[E] 



(A.6) 



l<a<fci 
Kb<N 



l<a<d 
Kb<N 



We define the sign of the imaginary parts of Kj\ by 
Im 

- Lj = J — . 

\imKji\ 

Integrating over E\, Eq. (1A.6j) reads 



(A.7) 



,. n (_ 1 \k 2 (k 2 -l)/2+(k 2 +l)k 1 N(N~l)/2 



x det 



no'-i)! 

3=1 



1 



Ber 



(A.8) 



K al — K fe2 J l<a<&i 
l<6<fc 2 



zL a / /(E^exphm^e^^l 

_/ I l<a<fci 



K6<JV 



{<2 1 } 



l<a<d 
l<b<k 2 



f(E)E 



6-1 







o-l 



| l<a<d 
1<6<JV 



In the lower right block we use the following property of the integral 
j f{E)E b - 1 (i-^ S(E)dE = for b > a . 



(A.9) 



Since d = N + k 2 — k\ < N, cf. Eq. (1131) . the last N — d columns in the lower right block 
in the determinant flA.8[) are zero. The matrix 



M 



/ ixm"- 1 (^) a s ( E ) dE 



(A.10) 



J Ka,b<d 
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is a lower triangular matrix with diagonal elements 

M ss = J f(E)E j ~ 1 ft AY 1 5 {E)dE = {-iy-\j - 1)! . (A.ll) 

R 

Thus, the determinant of this matrix is 

JV 

detM= (-^W-W (A.12) 

3=1 

We pull the Matrix M out the determinant (1A.8[) and find 



,_i^(*2-l)/2+(ij+l)fc. 



x det 
where 



l<a<fci 1 J l<a<fci 

l<6<fc 2 I H J d+l<b<7V 



i d r 

A" W («ai,«w) = V / /(E)E m - 1 exp[-m ol E]e(L S)dEM-i< 2 - 1 .(A.14) 

m,n=l K 

Again we use the fact that the determinant is skew-symmetric which allows also to write 



4/1 («) = ; = (A. 15) 



x det 



{^(«ai,«w)} 1<a<Jfci f(E)E b - 1 ex P [-tK al E]Q(L a E)dE 



l<a<fci 

l<b<fc 2 V R ) d+l<b<N 



for an arbitrary N G {d, d + 1, . . . , N}. For the cases (k 1 /k 2 ) = (1/1) and (ki/k 2 ) = 
(1/0), we identify 

Zy t («ol, ^62) = - K6 2 )A' (Ar) (K a i, K 62 ) , (A.16) 

= ~ lLa j fiE)^' 1 e M-^aiE]e(L a E)dE . (A.17) 

k 

This yields the result (fl4l) . 

Appendix B. Derivation of Eq. ( 1291) 

This derivation is similar to the one for the supermatrix Bessel function with Efetov- 
Wegner term in Sec. V.A of Ref. [TS]. We consider the integral 



gfflj (-l)"27r f r „. 3^ - V - J 

/"(0) (A-l)! 



^1/1 



e-^^-J e"^«J (r 2 ) d[p] . (B.l) 
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As in Ref. [18], we exchange the integration over the Grassmann variables by a 
differential operator which yields 

N-l 



(-1) 



N 



f(0) (N-l)i 

1 

1 



-tip 



+ 



d 



r\ — e l ^r 2 \dr\ 

The term 



+ e 



-up 



d_ 

dr 2 

d 



dr 2 

«i - K 2 /(ri) 



ri <9r 2 



«1 - K-2 
5 ; 

ri — e* v r 2 



exp[— «Str kt] 



(B.2) 



dr\dr 2 ■ 



(N - 1)! J ri - e^r 2 / (e^r 2 ) 



5 



exp[-zStr «r]rf ( e - **^- ) 5 (r 2 ) rfr lC /r 2 (B.3) 



N-l 



contains the supermatrix Bessel function with respect to U (1/1) [281 ESI HH] - The second 
term 



-1 



(iV-1)! 

i 

1 



rf e" 1 *-) 5(r 2 



d 



T\ — e l ^r 2 V dr\ 







+ e 



-iip 



dr 2 




N-l 



(B.4) 



dr- 



f (e^r 2 ; 



exp[— zStr nr] 



dr\dr 2 



T\ dr 2 

has to yield the Efetov-Wegner term. By partial integration, we evaluate the Dirac 
distribution and omit the generalized Wick-rotation. Thus, Eq. flB.2j) becomes 

'N-l 

(i\ — ij! / o> 1 " a j 1 * \ 

r\ 



x 



(iV-1)! 
f(n) 



lf(r 2 ) 



j=0 

exp[— tStr nr 



Qj+l Qj+1 

J- 1 \dridri dr J 2 +1 
dr\ . 



N 



,N-1 9 

1 «9rf 



(B.5) 



r 2 =0 



For all terms up to j — we perform a partial integration in ri and find a telescope 
sum. Hence, we have 

d 



dr\ 



fin) 
/(r 2 ) 



exp[— zStr «r] 



<iri = 1 



(B.6) 



r 2 =0 



This is indeed the Efetov-Wegner term. 
The second equality ( 129]) follows from 



I 1 (± + **±) 

J ri — e*^r 2 y dr\ dr 2 J 



N-l 



, i i /^rexp[-zStr «x]rf fe"^^- ) 5(r 2 ) 1 dndr 2 



.JV 



<9r 5 



iV-1 



ri — e*^r 2 \dri dr 2 J \f (e l ^r 2 ) 



exp[— zSti Kr] J dr\dr 2 



+ 



/ (e^r 2 ) 



exp[— «Str kt] 



N-l 



r\ — e 2 ^r 2 



d_ 

Or- 



N-l 



,N 



+ 



9r 5 



A? 



5(r 2 )dridr 2 
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r 



e -^__ ) 5 (r 2 ) 







+ (-1) 



N-l 



dr 2 

'N-l 



r\ — e % ^r 2 \dr\ dr 2 J \f (e l ^r 2 ) 



f(ri) 
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■ exp[— zStr nr] ] dr\dr 2 



N 



^ 1 dr{ f^ Q j\ 1 3r{ 



3=0 



fin) 

f (e^r 2 ) 



exp[— tStr nr] 



dri . 



r 2 =0 



Both sums cancel each other up to the term j = N which is the term l d N /dr^ in 
Eq. flR5j) . 



(B.7) 



Appendix C. Derivation of Eq. 

Let the characteristic function and, hence, the superfunction $( fcl / fc2 ) be factorizable, cf. 
Eq. (USD. To show the identity ([32]) we plug Eqs. (J2U]) and ((23} into the result ([HD for 
iV = d. We find 



Ber 



ft /(e*r i2 
Li=i 



-det d ri J ] | r 



-zi/> 



d 



Or 



J'2 



d-1 



<f(r i2 )d[r] 



(_l)(*:i+fc2)(fcl+fc2-l)/2^ 7r ^(fc 2 -A: 1 ) 2 /2-(fci+A : 2)/2 



2^, /Ber ^(/c) 



det 



{A(K 61 ,« a2 )} 



l<a<fc 2 
l<6<fci 



/(r 1 )rf- 1 e _t,Siiri dri 



d+l<a<Af 
l<6<fei J 



with 

A(K bl , K a2 ) 



exp(-2K W ri + ie^K b2 r 2 ] 
(«6i - «o2)(n - e^r 2 ) 



(C.l) 
(C.2) 



+ e" 



f(e^r 2 ) 







4 ( e"^— ) 5{r 2 ) 



d-1 



d[r] . 



The next step is to pull all factors of /, the monomials rf and the distribution 
(e~ l ^d j 'dr 2 j 5(r 2 ) out the determinant. Identifing the remaining terms of the 
integrands for all / we, then, get Eq. 



Appendix D. Derivation of the supermatrix Bessel function with 
Efetov-Wegner terms 



We use the result of Appendix C as an ansatz in Eq. (15U|) . To prove that this ansatz is 
indeed the result we are looking for we construct a boundary value problem in a weak 
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a 



\<a,b<k\ \<a<k\ 
l<b<k 2 

{Xab} Wab2} 

l<a<k 2 l<a,b<k 2 

K6<fei 



(D.l) 



with non-zero entries everywhere instead of a diagonal supermatrix k. Then the action 
of the differential operator 



Str 



da 2 



d 2 

a=l aaL l<a<6<fci 

+ 2 E 



daabidcr^ 



k 2 

E 

a=l 



5 2 



E 



a 2 



l<a<fe<fc 2 



d(Tab2d(T, 



,il,2 



l<a<k 2 
Kb<ki 



dx* ab dXab 



(D.2) 



on the left hand side of Eq. ( ]30l) yields 
<9 2 



Str 



da 2 



F(p) exp[— tStr ap]d[p] 



F(p)Str p 2 exp[— zStr ap]d[p] . 



(D.3) 



* 1 /fe 2 fcl/*2 

Since the integrand is rotation invariant the integral only depends on the eigenvalues of 
the supermatrix a. This leads to a differential equation in the diagonal supermatrix k. 
With the differential operator Sir d 2 /Ok 2 defined similar to Eq. ( PHI) we have 

I «2 



.(2) 
fci/fc 2 



d 

-Str— A/Ber^Li/.) 



Ber (kj 



(9k 2 



fci/fc 2 



exp[— «Str Kp]d[p] 



fc 1 /fc 2 



F(p)Str p 2 exp[— zStr Kp]c/[p] 



(D.4) 



;M0 



cf. Ref. EE]- 



The boundaries of Ej^L are given by E^ 1 ^ / fca _ 1 canonically embedded in Sj^L if 



one bosonic eigenvalue of a supermatrix in Ejj^L equals to a fermionic one, i.e. there 
are two numbers a 6 {1, . . . , fci} and 6 G {1, . . . , Z^} with K a \ = K^- For these cases 
we may use the Cauchy-like integral theorems for Hermitian supermatrices [T6l IMl fT8] . 
Without loss of generality we consider the and have 



F(p) exp[— zStr np\d[p] 



k\r\2—k\—k 2 



•l) fel 2 



F(p) exp 



-zStr k\ 



lK fci 1— K k 9 2 : 



g«0 

fe l/ fe 2 



-■fci-i/fca-i 



d[p). 

(D.5) 



Here we use the same symbol for the restriction of F on E^i 1 , fc2 _ 1 . 

The boundary condition (ID. 51) for the distribution (133|) can be readily checked. 
For the differential equation (1D.4f) we expand the determinant (I55j) in I < A; 2 rows and 
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columns in the upper block. Apart from a constant prefactor each term is given by 



Bex $(k,t, WBer[ 2) ,(r 



#(^ r ) = / ' ^ — 7^ II e ^ 6 ( r "iw) 5 ( r "2(m) 



J J exp(-iKo 1 ( )ir Wl ( )i) ]~J exp(ze lV ' ^(6)2^(6)2) , (D.6) 



x 

a=/+l a=Z+l 

where = diag (^(1)1, • • • , /Cwi(Qi, «wa(i)2, • • • , «w 2 (02) and r ^ = diag (r Wl(m)1 , . . . , r Wl(fcl )i, 
r uj 2 (i+i)2, ■ ■ ■ , r uj 2 (k 2 )i) with the permutations 0J\,U\ G &(kx) and U2,&2 G ©(^2)- The 
action of the distribution g(n,r) on Strr 2 is 



Ber^f^WBer^ ,(r. 



ff(«,r)Strr = [[e w d{r Wl[j)l )8{r W2{j)2 ) 

fel &2 

j [ exp(-?^ 1 ( a ) 1 r £Jl ( a ) 1 ) ] [ exp(ze # (6)2^(6)2) Strr 2 (D.7) 



fel &2 

X 

a=l+l a=l+l 

because all other terms are zero due to the Dirac-distributions. The differential operator 
in Eq. (1D.4|) acts on g(n,r) as 

1 Str ^/ier^ 

Ber£ 2 > (k) 



Str ^V Ber *iW K ^ (K ' r) 



52 WBerf^Kflj) A/Ber^_ z/ (r w ) ' 



Berg fe2 (K) Ber^(r) 

fel &2 

x I [ exp(-z^ 1 ( a) ir Wl ( a) i) J J exp(?e # ^2(6)2^2(6)2) ■ (D.8) 

a=l+l a=l+l 

We split the differential operator Strd 2 /dK 2 into a part acting on kq and a part for 
the remaining variables in k. The term for the latter variables acts on the exponential 
functions in Eq. flD.8j) and contributes the term —Strr 2 . For the term according to k$ 
we use the identity 

Str J^Ber^HO. (D.9) 

Thus, the differential equation is also fulfilled by (ph x /k 2 - 

Appendix E. Double Fourier-transform 

We consider the integral 

-2Tre~ l ^S(r b i)S(r a2 ) exp (-in bl r b i + in a2 r a2 ) , _^ 

H — h X\^b\ - e ^K a2 



det 



Kn - e % ^K a2 r bl - e^r a2 J i<a<k 2 

l<b<fcl 

{r^exp (-iK hl r bl )\ 

l<6<fei 
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x det 



f 2ire % ^5(s bl )5(s a2 ) _ exp (ir bl s bl - ir a2 s a2 ) 



r bl - e^r a2 



Sbi ~ e-»^So2 



X(ni - e # r a2 ) 



Ki 1 exp(«r 6 is 61 )} 



l<a<fc 2 
l<fe<fci 



k 2 



l<fe<fci 



X 



d[r] 



We omit the two sums over the permutation groups, see Eq. f l3"3"j) . They do not contribute 
any additional new information of the calculation and the missing terms can be regained 
by permuting the indices of the eigenvalues in s or k. 
The expansion in the first determinant yields 

k 2 



(E.l) 



J = E £ 



sign UiU 2 



1=0 wiee(ki) 
^ 2 e6(fc 2 ) 



(Z!)2(fc! -Z)!(A: 2 -Z)! 



det 



l+*2 



-2?re *^'^(r tJl ( i ,)i)J(r CJ2 ( a )2) 



x det 



exp (-«« £Jx (6)l'"a ) x(6)l + WC W a(o)2r a « 1 (a)2) , _w, s 



{C(l)i eX P H^i(fc)l r -i(6)l 

f 27re"^(s b i)£(.s Q2 ) _ exp (zr b is b i - ir a2 s a2 ) 



x det 



^ } l<a<fci 



J l<a,b<l 



l+l<a<k 2 
2+l<6<fci 



k 2 

l+l<b<kx 

lip 



r b i - e^r a2 



s bl - e l ^s a2 



X(r b i - e^r a2 ) 



{s a bl x exp (ir b is bl )} 



l<a<k 2 
l<fe<fci 



l<a<fci— k 2 
l<6<fci 



X 



d[r] 



sign uiu 2 



x 



b <W^V b </^ b <a> 

(_l)fci(fcl-l)/2^ 

w 2 G6(fc 2 ) 

fci k 2 

-tip 



(unh-mh-iy. 



det 



-(27r) 2 e^^(^ l(fc)1 )^(^ 2(a ) 2 ) 



J Ko,6<! 



fci fc 2 „ fcl fe2 

JJ JJ x(«wi(6)i - e" iV, ^ 2 ( a ) 2 ) / JJ exp (-m Wl(a) ir a i) JJ exp (iK^^ar^) 

a=Z+lft=Z+l „i.,xl„-9; a=^+l b=l+l 



x det 



2v:e l " p 5{s UJl ( b )i)5{s L0 ^ a)2 ) exp ^r^s^^i - zr a2 s W2(a)2 ) 



r bl - e^r a2 



S Wl (6)l - e l ^S W2 (a)2 



X(r u - e"V o2 ) 



Z+l<a<fc 2 
/+l<6<fci 



fc 2 

/+l<6<fei 



X 



d[r 



Beri 2) „ (relBerPl, (si 



(E.2) 



ki/k 2 \ I ki/k 2 y 

where the function "sign" yields 1 for an even permutation and —1 for an odd one. The 



On the Efetov-Wegner terms by diagonalizing a Hermitian supermatrix 



24 



permutations in the indices of the r are absorbed in the integration. We remark that 
the remaining integral goes over ki + k 2 — 21 variables because we have already used the 
Dirac-distributions. 

With help of the formula 



2ire li '5(s UJl{b)1 )5(s 

LU2(a)2) 

r&i - e l ^r a2 

2me t *5(8 Ul (fi )1 )5(8 u] 2 (a)2 ) 
K 2 («wi(6)i - e-^K u ^ a p)(r bl - e^r a2 ) 



exp ( )2r a 2) d[r] 



d 



+ e 



-zip 



d 



dr bl dr a2 
(2n) 2 e^5( Su)im )5( Su)2ia)2 ) 



exp (-iK Ul(b)1 r bl + m W2{a)2 r a2 ) d[r] 



we integrate and sum the expression (1E.2j) up. This yields 



^fcx/fca 0«, r)fi kl/k2 (-ir, «)Ber g } /fca (r)d[r 



(E.3) 



(E.4) 



k 1 + k 2 



/_ 1 \fci(fci-l)/2 7r (fo-fci) 2 



2 2fc! fc 2 -fa -fa fcl | k2 1 ^ / Ber (2)^ , 



«0 Ber fa 



E det 



/A . 2 (sju, ie e(fa) 

w 2 e6(fa) 



a 6 



^1(5)1^(^61 - 



Mb)l; 



with 



a6 



e # 5(s a , 1 ( 6 ) 1 )5(s a)2 ( a)2 ) 



«61 



^a2 



[1 - X(«61 



Ka2)] 



+ e —3; Xi^bi - e r /? a2y 



(E.5) 



This is the result (1411) . The index a goes from 1 to k 2 in the upper block and from 1 to 
k\ — k 2 in the lower block whereas b takes the values from 1 to k\ in both blocks. 

Appendix F. Calculations for subsection 16.11 



We diagonalize the supermatrices a and p in Eq. (T49i) and have for the generating 
function 

1 f d[r]d[s] 



4 N v 



(2ni) 2k 
x det 



:J-$( S )/ ; 



(F.l) 



2vre ^(r^Wr^) exp {-iK bl r bl + le^K^r^) 

H ~ - K a2j 



r 61 - e^r a2 



x det 



K b l — K-a2 

2'Ke^5(s b i)5(s a2 ) exp {ir bl s bx - w a2 s a2 ) 



Ka,6<fc 



r 6 i - e^r a2 



Sbl 



-lip 



Sa2 



X(r bl - e^r a2 ) 



l<a,b<k 



This expression is not well defined because the supersymmetric Ingham-Siegel is at zero 
not well defined. We recall that the supersymmetric Ingham-Siegel integral factorizes in 
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each eigenvalue of the supermatrix r, cf. Eq. f )50|) . To understand Eq. (IF.lj) we have to 
know what (0) is. Since the supersymmetric Ingham-Siegel integral is a distribution 
we consider an arbitrary rotation invariant, sufficiently integrable superfunction / on 
the set of (1 + 1) x (1 + 1) Hermitian supermatrices. Then we have 



/ 



/ 



f{ P )l[ N \p)d[p] 



1/1 V+/1 



f(p) exp(— iStr pa + eStr p)d[p] 



Sdet ~ N (a + tet 1+1 )d[a] 



J 



1 / /(p)exp(£Strp)d[p] 



= /(0) 

= -z/^/f^O) (F.2) 

with help of the Cauchy-like integral theorem for (1+1) X (1+1) Hermitian supermatrices, 
see Ref. [16, 18J. Please notice that the constant resulting from the Cauchy-like integral 
theorem converts to the complex conjugate when the generalized Wick-rotation is 
complex conjugated. The last equality in Eq. (1F.2|) is the Cauchy-like integral theorem 
formally applied to the left hand side of Eq. (1F.2[) . Hence we conclude that 1^(0) = % 
in a distributional sense. Using this result we find 



x det 

x det 
x det 



(_l)fc(fc+l)/2 k 



2k(k+l) 7T 2k 



E E 

1=0 0Ji,w 2 GS(fc) 



sign UJ1UJ2 



(F.3) 



exp (-m Wl(6)1 r Wl(fc)1 + ze 4,/, /t aJ2(a)2 r W2(a)2 ) {JV) 

=7 H +^(6)1, e V r Wa(o ) 2 )x(«wi(6)l - «wa(o)2j 



27re # (5(s fe i)5(s a2 ) exp (zr w s w - ir a2 s a2 ) 



J Z+l<a,fc<fc 



r 6i - e^r a2 



Sbl 



- P -%il> 



Sa2 



-1 



-27ae # ^(r a;i ( & )i)(5(r a;2 ( a)2 ) 



x( r w - e^r a2 
d[r]d[s] 



J Ka,b<fe 



i<a,fe</ WBer 



'SJ r )v Ber ^ 



signa;ia; 2 



O/.M/.-- I) 7r 2fe' zJ /|.;. II 

Z=0 wi,u) a e6(fc) L 1 -^ V-J a ,6=i+l 



J ) X( K a)i(6)l - «wa(a)2) / 



R 4fe-2I 



X 



x det 



x det 



] J exp {-iKu^mi + ^e # K tt , 2(a)2 r a2 ) ii (r w , e # r a2 ) 
1 

27re^5(s Wl(6 )i)5(s W2(a ) 2 ) exp (irus^^i - «r a2 s a , 2(a)2 ) 



a,6=Z+l 



r H - e^r a2 



s Wl (6)i - e # s W2 ( a ) 2 



- e # r a2 ) 



Z+l<a,fe<fc 



-(2?r) 2 ze^^(5 a;i (j ) )i)^(5 a;2(a) 2) 



d[r]d[s] 



J / v /Ber^ fc (/c) 
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1 k 



sign U1OJ2 



x 



2fc(fc— i)^ 

d[r]d[s 



S(fc) o,6=i+l 



I [ X(«wi(6)l ~ (a 



,(2) / 
fc/fc 1 ' 



,J*$(s) det 



Ber ):;,.{ k) 



x det 



e' V '^(So J1 (6)l)^(s^( a )2) 
^6Ji(fe)l K w 2 (a)2 



■Sa; 2 (a)2j 

K aj 2 (a)2 



)2j 



Ka,b<l 



N (e # s W2 (a)2 - K U2{a)2 ) 



N-l 



2m (e ^s Wl(6 )i - « Wl ( 6 )i) JV+1 - e »^s^( a ) 2 ) J 

We perform the sum and use the identity 



K 



w 2 (a)2 



K 



N 



X( K LOi(b)l - ^2(0)2) 



cj 2 (a)2 



Wl (6)l 



A' 



(F.4) 



Then we have the result (EI 



Appendix G. Calculations for subsection 16.21 

In the first step we derive the Fourier-transform of the superdeterminant in Eq. (|52|) . 
Let the entries of the diagonal (k + k) x (k+k) supermatrix r and the entries of diagonal 
N x N matrix Eq be the eigenvalues of the supermatrix p and the Hermitian matrix 
H , i.e. p = UrW with U e V (k/k) and if = V^V^ with V G U(iV). Then the 
Fourier-transform is 



J= y exp(-zStrrcr + eStrp)Sdet l {a ® 1 N + at k+ k ® H + let 
'2ne^S(s bl )S(s a2 ) exp (-is bl r bl + is a2 r a2 



-'k/k 



2 fc2 7T fc 



det 



N(k+k)j 



X(ni - e # r a2 ) 



r 61 - e^r a2 



Sbi - e-^So2 



J Ka,b<k 



x Sdet x (s ®t N + alfe+fc ® # + iet N {k+k)) 



d[s] 



(G.l) 



With help of identity (jSJ) we find 
z fc exp(eStrr) ^-^ 



sign U1U2 



x 



x det 



J Sdet x (s ® Ijv + al k+ k ® E + tel N(k+k)) det 
exp (-zs bl r Ul (p)i + 2s a2 r W2 ( a ) 2 ) 



2Tre^S(s bl )S(s a2 ) 



d[s] 



Z+l<a,fe<fc 
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exp(eStr r) 



2 k2 n k 



x 



det 



E E 

=0 £ji,0)26S(fe) 
k 



sign U\U 2 



A N (aE )^BerfJ k (r) t^o Ul ,^e@(k) ^ k ^ 0,6=/+ 1 

2ne^5(s bl )5(s a2 ) 
r Ul (b)i - e*^r W2 ( ) 2 



II exp (^^tfja) 



(G.2) 



l<a,b<Z j = j_(_l 



x det 



exp (-zs b ir Wl ( fe )i) 




-^s a2 + ie\ 


N 


{e-^ + ze)"- 1 




V Sbi + 1£ ) 


exp (-2Sbir Wl(fe) i) 




N 


(-aE^r 1 


Sbi + «e + a-Ea 0) 


ys bl +ie J 



d[s] 



In the left upper block both indices a and b run from I + 1 to k whereas in the right 
lower block the range is from 1 to N. In the right upper block a goes from I + 1 to k 
and b goes from 1 to N whereas it is vice versa in the left lower block. We sum all terms 
in Eq. (I0.2[) up and pull the integrations into the determinant. Then we have 



J 



1 



2 k2 n k 



A N (aE ) x /Bev k %(r 



det 



{Mn^r a2 )}^ ah< _ k {A b2 (e^r a2 )} ^ 

l<b- 

Wn.,««l, {(-"WU. 



Kb<N 



l<a<N 
Kb<k 



b<N 



(G.3) 



where 

v4i(rbi,e # r a2 ) = exp [e(r M - e # r a2 )] 



2tt 



r 6 i - e^r a2 



+ 



exp {—is x r bl + is 2 r a2 ) ( e ^ s 2 + is 
S\ — e~ l ^s 2 \ si + is 



A' 



X(ni - e l1p r a2 )d[s 



- 2ni 



l[ N) (r bl ,e^r a2 ) 
r bl - e^r a2 
{2-nfi <8(r 61 ) 
(N -l)\r bl -e^r a2 



-itp 







Or- 



al 



JV-1 



e"^5(r a2 ) 



(G.4) 



Ab 2 (e^r a2 ) 



exp (-£e # r a2 ) / exp (zs 2 r a2 ) (e~^s 2 + isf 1 e~^ds 2 J]^i - e # r a2 ) 



27r ( — ze" 



5 



b-l 



A 3 (rbi,a^ 0) ) = exp(erbi) 



dr a2 

exp ^-zsir 6 i 



i=i 



(G.5) 



N 



s 1 + is + aE a 0) \si+ie 
,(0) >„ fc 

a ' 

n! 



n=7V 



dsiflx(r6i-e^r i2 ) 



j'=i 



(G.6) 
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Surprisingly this part of our result agrees with the one in Ref. [10J (apart from a forgotten 
2txi in the upper left block of Eq. (6.9) in Ref. [TU] and the characteristic functions 
xi r bi — e"V a 2)) although we omitted all Efetov-Wegner terms in this work. 
The second step is to diagonalize the supermatrix p in Eq. 



Z { k N \K,aH ) 



{ _ z) k(k+i)/2 / det 



(G.7) 



-2ire ^5(r bl )5(r a2 ) exp (-iKbim + «e # « a2 r o2 ) , s 
H X(«6i - K a2 ) 



Kbl — Ka2 



r bl - e^r a2 



Ka,b<k 



x det 



l[ N \r bll e^r a2 ) 






6-1 k 

e-^5(r a2 )l[x(r n -e^r a2 ) 

i=i 


2n(r bl - e^r a2 ) 


\ 6 dr a2 ) 


^ {iaE { a ] r bl ) n # 
e ( r w)2_, n \ [[x(ni e^r j2 ) 

n=N ' j=l 





$(r)d[r] 



A iV (aE );/Beri 2 / ) fc ( fi :)JBer( 2 / ) fc (r) 



The range of the indices in the second determinant is the same as in Eq. flG.3j) . 
Expanding the first determinant we have 



zl N) (K,aH ) 



(Gi 



_ % \k{k+X)/2 



exp {-iKu^r^h)! +ie l ^K LLj2{a)2 r ul2{a)2 ) 



E E 



sign wia;2 
[*!(fc-0!]\. 



x 



]Q x( K a;i(6)l - «w 2 (a)2) 



det 



J Z+l<a,6<fc 



x det 



x det 









6-1 fc 

e- 4 ^(r a2 )n^(^i- e ^) 
i=i 


27r(r 61 - e^r a2 ) 




e (rH )E ( ;, m) 

n=iV ' j=l 





-27re- l ^(^ l(6)1 )5( 



)2j 



1^1(6)1 — K aj 2 (a)2 



$(r)d[r 



1 / jBerj^r 



When integrating and summing up we use the normalization (IF. 2D of the supersymmetric 
Ingham-Siegel integral and arrive at the result fl53|) . 
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